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Abstract 

Using an effective interaction approach to describe the interactions between the dark matter par- 
ticle and the light degrees of freedom of the standard model, we calculate the antiproton flux due to 
the annihilation of the dark matter in the Galactic Halo and compare to the most recent antiproton 
spectrum of the PAMELA experiment. We obtain useful constraints on the size of the effective 
interactions that are comparable to those deduced from collider and gamma-ray experiments. 
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I. INTRODUCTION 



The presence of cold dark matter (CDM) in our Universe is now well established by a 
number of observational experiments, especially the very precise measurement of the cosmic 
microwave background radiation in the Wilkinson Microwave Anisotropy Probe (WMAP) 
experiment [1]. The measured value of the CDM relic density is 

tt CDM h 2 = 0.1099 ±0.0062 , 

where h is the Hubble constant in units of 100 km/Mpc/s. Though the gravitation nature 
of the dark matter is established, we know almost nothing about its particle nature, except 
that it is nonbaryonic and to a high extent electrically neutral. 

One of the most appealing and natural CDM particle candidates is weakly-interacting 
massive particle (WIMP). It is a coincidence that if the dark matter (DM) is thermally 
produced in the early Universe, the required annihilation cross section is right at the order 
of weak interaction. The relation between the relic density and the thermal annihilation 
cross section can be given by the following simple formula [2] 

Olpb 

x ~ (av) ' 1 j 

where (av) is the annihilation rate of the dark matter around the time of freeze-out. Given 
the measured f2cDM^ 2 the annihilation rate is about 1 pb or 10~ 26 cm 3 s _1 . This is exactly 
the size of the cross section that one expects from a weak interaction process and that would 
also give a large to moderate production rate of the WIMP at the Large Hadron Collider 
(LHC). In general, production of dark matter at the LHC would give rise to a large missing 
energy. Thus, the anticipated signature in the final state is high-pr jets or leptons plus 
a large missing energy. Note that there could be nonthermal sources for the dark matter, 
such as decay from exotic relics like moduli fields, cosmic strings, etc. In such cases, the 
annihilation rate in Eq. ([Q) can be larger than the value quoted above. 

There have been many proposed candidates for the dark matter. Instead of specifying a 
particular model we adopt an effective interaction approach to describe the interactions of the 
dark matter particle with the standard model (SM) particles 0-8]. One simple realization is 
that the dark matter particle exists in a hidden sector, which communicates to the SM sector 
via a heavy degree of freedom in the connector sector. At energy scale well below this heavy 
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mediator the interactions can be conveniently described by a set of effective interactions. 
The strength of each interaction depends on the nature of the dark matter particle and the 
mediator. An interesting set of interactions between the fermionic dark matter \ and the 
light quarks q can be described by (qTq)(x^'x)i where r,T' = a^ v , a^'y 5 , 7^, 7 M 7 5 , 7 5 and 
1. Note that due to the following identity 

the axial tensor cr^7 5 is related to the tensor a a/3 and thus should not be regarded as an 
independent set. A more complete set of interactions involving fermionic and scalar dark 
matter candidates that we will study in this work are listed in Table [H Without a particular 
model in mind we will treat each interaction independently in our analysis. 

There have been some recent works on constraining the interactions at present and future 
collider experiments {3- 5] and using gamma-ray experiments Fan et al. {t| also wrote 
down the effective nonrelativistic interactions between the dark matter and nuclei. There 
was another work in which the dark matter couples only to the top quark and corresponding 
predictions at direct and indirect detection experiments as well as colliders are obtained js]. 

In the present work, we first estimate the lower bounds on the new interactions based on 
the fact that if a particular interaction is the only contribution that can thermalize the dark 
matter particle in the early Universe, we require this interaction must be strong enough such 
that the resulting relic density would not overdose the Universe. In addition, we proceed 
to calculate the antiproton flux coming from the effective interactions. We expect the latest 
antiproton data from PAMELA [9( can put a strong constraint on the size of the interactions, 
based on the fact that the existing data do not allow excessive flux above the conventional 



background. On the other hand, the positron spectrum from PAMELA 



10j showed some 



excessive above the conventional background and thus if we used it to constrain the model, it 



in general gives a weaker constraint than the antiproton flux [8l 111]. Therefore, in this work 
we focus on the constraints from antiproton flux. Similar ideas of using the antiproton flux 
from earlier PAMELA data was considered in Ref. 12] to confront the low energy CoGeNT 
experiment [13] for a lower mass DM. 

The organization of the paper is as follows. In the next section, we describe the inter- 
actions between the dark matter particle and the SM particles, in particular quarks and 
gluons. In Sec. Ill, we study the velocity dependence of the effective operators based on 
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the nonrelativistic reduction. In Sec. IV, we calculate the annihilation rates during the 
freeze-out and make sure that the interactions would not overdose the Universe. In Sec. V, 
we calculate the antiproton spectrum due to the dark matter annihilation in Galactic Halo. 
We compare with other constraints and conclude in Sec. VI. 



II. EFFECTIVE INTERACTIONS 



Let us start by assuming the dark matter is a Dirac fermion and its effective interactions 
with light quarks via a (axial) vector-boson or tensor-type exchange are given by the following 
dimension 6 operators 

C 

A=i-6 = -r| (xTiX) (<?r 2 g) , (2) 

i 

where r 12 = 7 M , 7^7 5 , o~^ u or a Ml/ 7 5 with a^ v = i{pi llr f — 7^7 M )/2, and C is an effective 
coupling constant of order 0(1). For Majorana fermion the I\ = 7^ or type interaction 
is identically zero, and so for vector or tensor type interaction the fermion x i* 1 Eq.© is 
understood to be Dirac. Explicitly, we assume the dark matter candidate to be Dirac, but 
the results are also applicable to Majorana dark matter. 

Next set of operators are associated with (pseudo) scalar-boson-type exchange 

£i= 7 -io = ^tit (x r ix) , (3) 

i 

where r 12 = 1 or 27 s . The m q dependence in the coupling strength is explicitly shown for 



scalar-type interactions. We use the current quark masses in the Lagrangian given by 

m u = 0.0025 GeV, m d = 0.005 GeV, m s = 0.101 GeV, 
m c = 1.27 GeV, m h = 4.19 GeV, m t = 172 GeV. 

Another light degree of freedom that couples to the Dirac dark matter is the gluon field 1 

4=n-i a = Ca f K T x) (XT X ) G^G« U (4) 

i 

A=i3-i4 = Ca f™ x) (XT X ) G^G% (5) 



1 We do not study the other gauge bosons, like W and Z bosons, because they decay into light quarks 
which then fragment into p. The secondary antiproton spectrum would be softer in this case. 
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where T = 1 or 27 s and the strong coupling constant is evaluated at the scale 2m x where 
m x is the dark matter mass. 

Finally, we also write down the corresponding operators for complex scalar dark matter. 
Again, we note that the interactions for real scalar dark matter is similar to complex one 
and differ by a factor of two. We simply focus on the complex scalar dark matter. The 
operators corresponding to vector boson exchange are 

£=15,16 = ^ (x^Xx) (qi^q) , (6) 

^ y 

where T = 1 or 7 5 and x^o^X = xK®nX) ~ (^uXOx- Those corresponding to a scalar boson 
exchange are 

A=ir,i8 = ^ (x f x) (qVq) , (7) 

i 

where T = 1 or ry 5 . The corresponding gluonic operators are 

A=i9 = ° a f™ x) (X f x) G^G a , v , (8) 

i 

= lCa f ^ ( x t x ) G ^Q ativ . (9) 

The whole list of operators are listed in Table [I] We will consider one operator at a time, 
and set the coefficient C = 1 for simplicity. 

Note that in calculating the annihilation rate in the freeze-out in the early Universe, we 
include all light-quark flavors (u,d,s,c, b) as well as the heavy top quark which is relevant 
when m x rises above the top quark threshold. However, in the calculation of the antiproton 
flux from dark matter annihilation in the present Galactic halo, we only include the light- 
quark flavors. We ignore the \X ~^ ti contribution, because the t and i first decay into 
bW — > bqq' before each light quark undergoes fragmentation into hadrons, including proton 
and antiproton. Therefore, the antiproton spectrum would be significantly softer than the 
direct fragmentation as in \X ~ * |S|] . We anticipate that by ignoring the ti contribution 
the limits we obtain from the PAMELA data would not be affected to any significant amount 
in the case we just use the five light-quark flavors. 

III. VELOCITY DEPENDENCE IN THE NONRELATIVISTIC LIMITS 



In order to easily understand the results that we obtain in Sec. V, we are going to examine 
the dependence of the annihilation cross section on the velocity of the dark matter particle 
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TABLE I: The list of effective interactions between the dark matter and the light degrees of 
freedom (quark or gluon). We have suppressed the color index on the quark and gluon fields. 
These operators have also been analyzed in Refs. (^,0,13]. 



Operator 


Coefficient 


Dirac DM, Vector Boson Exchan; 


ge 


01 = (x7 M x) (qi^q) 

02 = (X7 M 7 5 X) (QIixQ) 

03 = (x7 M x) (<fy M 7 5 <?) 

4 = (X7 M 7 5 X) (<77^7 5 «) 

5 = (x^x) (q^ug) 

Oe = (x^7 5 x) (?<V<z) 


A 2 " 

c 

A 2 " 

C 

a 2 " 

c 

A 7 

o 
A 7 

C 
A 2 


Dirac DM, Scalar Boson Exchanj 




07 = Odd (qq) 

08 = (x7 5 x) (qq) 

09 = (xx) (qi 5 q) 

010 = (x7 5 x) (qi 5 q) 


Cm q 
A3 

A^ 

A3 
Cm, 
A« 


Dirac DM, Gluonic 


On = (xx)G^ 
Oi2 = (x7 5 x) G^G^ 


4A 3 " 

iCa s 
4A» 



013 = (XA)0^0^ & 

014 = (X7 5 X)0^0^ gfc 

Complex Scalar DM, Vector Boson Exchange 

015 = (X f 5x) (97^) § 

Pie = (x ] Kx) (yrVg) %_ 

Complex Scalar DM, Scalar Vector Boson Exchange 

On = (X f x) (qq) ^ 

Oi8 = (X f x) (qiS) 

Complex Scalar DM, Gluonic 



019 = (X f x) G^W 

020 = (X f x) G^G^ 



Ca s 
4A 2 " 

iCa s 
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in the nonrelativistic limit. The current velocity of the dark matter in the Universe around 
the Sun is about v 300 km s -1 10 _3 c, where c is the speed of the light. Given such a 
small v the dependence on v is very important. For instance, the annihilation rate of the 
Dirac DM with a scalar boson exchange, given by the interaction in O7, would suffer from a 
factor of v 2 . Therefore, we expect the antiproton flux from such an operator would be very 
small. 

Let us consider the operators 0\ to 6 with (axial) vector-boson/tensor-like exchange. 
In terms of Dirac spinors (if) and ip) the relevant part of the annihilation amplitude of the 
Dirac DM is given by 

<Mp 2 )r>(pi) (10) 

for x(P2)x(Pi) ~~ ^ QQ an d T = 7^, 7^7 5 , a^ u or a-^7 5 . In Dirac representation, the gamma 
matrices are given by 

'o / 
/ 

where <7j(i = 1,2,3) are the Pauli matrices. In the nonrelativistic limit, the spinor for the 
DM x is ^ — £ ( ) ? where e — 0(v/c). On the other hand, the spinor for the antiparticle 

w 

X is if) = rf(e, 1)7°. 2 Therefore, we can expand ip^ip, in the nonrelativistic limit, as 

■07°'0 ~ 2e?7'^ 

^7> ~ (1 + e 2 )^^ 

where the space-like parts are not suppressed by v/c. On the other hand, ip^^ 5 ip in the 
nonrelativistic limit are 

^7°7 5 V> ~ (1 + e 2 )^ 
■?/>7*7 5 -?/> — 1etfo£ 

where the space-like parts are suppressed by v/c. It is clear that in the nonrelativistic 
limit the time-like and space-like parts behave very differently. We can then consider them 
separately when it is squared, traced, and contracted with the trace of the light quark leg. If 



It is different from the direct scattering with a nucleon, where we need the = £^(1, e)7°- 
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we look at the trace of the part {q^y^q) or (q^^q) in the annihilation amplitude, the time- 
like part after being squared and traced gives a quantity close to zero, while the space-like 
part after squared and traced gives a quantity in the order of vr? x . Therefore, it is clear now 
that ^7^ multiplied to (q'J^q) or (qj^q) will not be suppressed, while i/jj^^ifj multiplied 
to {q^f^q) or (qj^q) will be suppressed. The above observation is consistent with the results 
that we obtain in Sec. V. From Table [Til the limits on 0\ and O3 are much stronger than 
those on 2 and O4. The operators O5 and Oq contain unsuppressed components in fiu = Oi 
entries. 

In contrast, the operators 7 to O w with (pseudo) scalar-boson exchange are suppressed 
when there is no 7 5 in the fermion line of x> which is obvious from the following in the 
nonrelativistic limit 

■07 5 ^ ~ . 

Again, it is then obvious from Table [TTJ that the limits on Og and 0\q are much stronger 
than those on O7 and Og. The gluonic operators in On_i4 follow similar patterns: 011,13 
are suppressed relative to Ou,u- 

It is also straightforward to understand the velocity dependence for the scalar DM, rep- 
resented by the operators Oi5_2o- Except for Oi^iq all of them are not suppressed by v/c, 
because of the presence of the <9 M in 015,16. This will bring down in the vertex factor. 
While p Q is of order m x , is v/c. Therefore, when it contracts with the quark leg, the 
overall result is suppressed by v/c. 

IV. ANNIHILATION RATES AROUND THE FREEZE-OUT 

It is obvious from Eq. ([T|) that if the annihilation rate falls below 1 pb, then the thermal 
relic density would be more than the WMAP data can allow. Therefore, we have to restrict 
the annihilation rate to be larger than about 1 pb. More precisely, using the most recent 
WMAP result on dark matter density ficDM^ 2 = 0.1099 ± 0.0062 [l[ the annihilation rate is 

(av) ~ 0.91 pb . (11) 

We assume v ~ 0.3 at around the freeze-out time in the early Universe. 
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We calculate the annihilation rates for all the operators and show the contours in Fig. [TJ 
Fig. [2j and Fig. [3] for Dirac DM with (axial) vector-boson/tensor-like exchanges, Dirac DM 
with (pseudo) scalar-boson exchanges and Dirac DM with gluonic interactions, and scalar 
DM, respectively. The solid lines are the contours in (m x ,A) plane with annihilation rate 
av = 0.91 pb. All the values of A above the solid lines would give a too small annihilation 
rate, and thus would result in a too large thermal relic density. Therefore, the A below the 
solid lines is the allowed region. The cusp structures in the plots are due to the onset of the 
top quark contributions when m x > m t . The dashed lines are the limits from the antiproton 
flux, which will be explained in the next section. 

6000 
5000 
4000 

> 

C33000 

< 

2000 



1000 



100 200 300 400 500 

m x (GeV) 

FIG. 1: The upper limits on Aj due to the relic density constraint, requiring av > 0.91 pb for 
operators Oi_6 involving Dirac DM with (axial) vector-boson/tensor-like exchanges (shown by 
solid lines). The lower limits on Aj due to the antiproton- flux constraint at 3a level for the same 
operators (shown by dashed lines with the corresponding color). 
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FIG. 2: The same as Fig.[TJ but for operators O7-14 involving Dirac DM with (pseudo) scalar-boson 
exchanges (O7-10) and Dirac DM with gluonic interactions (On-u). 

V. ANTIPROTON FLUX 

An important method to detect the dark matter is by measuring its annihilation products 
in Galactic halo. Current experiments can detect the positron, antiproton, gamma ray, and 
deuterium from dark matter annihilation. The Milky Way halo may contain clumps of dark 
matter, from where the annihilation of dark matter particles may give rise to large enough 
signals, such as positron and antiproton, that can be identified by a number of antimatter 
search experiments. The most recent ones come from PAMELA [9|, I10J, which showed a 
spectacular rise in the positron spectrum but an expected spectrum for antiproton. It may 
be due to nearby pulsars or dark matter annihilation or decays. If it is really due to dark 
matter annihilation, the dark matter would have very strange properties, because it only 
gives positrons in the final products but not antiproton. Here we adopt a conservative 
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FIG. 3: The same as Fig. [TJ but for operators O15-20 involving complex scalar DM. 



approach. We use the observed antiproton as a constraint on the annihilation products in 
XX annihilation. In general, the positron data would give a weaker constraint as it allows 
some level of signals of dark matter annihilation jg]. 

The antiproton flux can be obtained by solving the diffusion equation with the correspond- 
ing diffusion terms and the appropriate source term for the input antiproton spectrum: 

2 



Q, 



v 



PCDM 

Mr 



W \ dN P 



(12) 



'CDM / 

where 77 = 1/2 (1/4) for (non-)identical initial state, and Tp is the kinetic energy of the 
antiproton which is conventionally used instead of the total energy. We again solve the 
diffusion equation using GALPROP 15] . 

In our case, the dominant contribution comes from 



xx -> qq -> v + x 



(13) 



11 




FIG. 4: Antiproton fraction spectrum predicted for the operator 0\ = xy(x7^x) (<?7m?) f° r a f ew 
values of A. The mass of the dark matter is chosen to be 200 GeV here. The data points are from 
PAMELA 

in which all the q, q (q = u, d, c, s, b) have probabilities fragmenting into p. We adopt a 
publicly available code [l6] to calculate the fragmentation function D q ^h{z) for any quark 
q into hadrons h, e.g., p,p,ir. The fragmentation function is then convoluted with energy 
spectrum dN/dT of the light quark to obtain the energy spectrum of the antiproton dN/dTp. 
The source term dN/ dTp is then implemented into GALPROP to calculate the propagation 
from the halo to the Earth. We display the energy spectrum for the antiproton fraction in 
Fig. H] for the operator Ox, in which various values of A are chosen. We only chose a typical 
operator. The effects of other operators are similar. 

Here we adopt a simple statistical measure to quantify the effect of each operator. We 
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calculate the 3a limit on each scale Aj. We assume the data agree well with the expected 
background, and then we calculate the \ 2 with finite Aj's until we obtain a \ 2 difference 
°f X 2 ~ Xbkgd = 9 (3c)- Note that the uncertainties in the background estimation of the 
low energy range (< 4 GeV) are large in GALPROP, mainly because of different profiles 
employed. We therefore focus on the data points above 4 GeV when we calculate the \ 2 ■ 
The data points above 4 GeV enjoy a small \ 2 — 5-0 for 13 degrees of freedom. We tabulate 
all the lower limits of A^s in Table U for m x = 50, 100, 200 and 400 GeV. 

Limits for intermediate values of m x are shown in Fig. [TJ Fig. |2j and Fig. [3] for Dirac 
DM with (axial) vector-boson/tensor-like exchanges, Dirac DM with (pseudo) scalar-boson 
exchanges and Dirac DM with gluonic interactions, and scalar DM, respectively. The solid 
lines are the upper limits due to thermal relic density. The dashed lines are the lower limits 
due to antiproton flux. Therefore, for each operator there is a valid range of A. For example, 
the operator Oi >3 requires 1.6 TeV < Ai j3 < 3 TeV for m x = 200 GeV. The best limit is 
obtained for the Dirac DM with tensor interactions in 0$^ at 1.9 TeV < As^ < 3.6 TeV for 
m x = 200 GeV. In general, the Dirac DM with vector-like exchanges gives the best limit, 
except for the operators 02,4, which are well known for velocity suppression. The operators 
Oj^io for Dirac DM with scalar-boson exchanges naturally give rather weak limits, because 
of the factor m q in the coupling constant. In addition, operators 07,9 are suppressed by the 
velocity. The gluonic interactions in the operators 0n_i4 also give mild limits because of the 
a s ~ 10 _1 in the coupling constant, in which the operators On,i3 are further suppressed by 
velocity. On the other hand, the operators for scalar DM give rather mild limits, especially 
015,16 give the weakest limits, because the derivative couplings in Oi5,i6 bring down a factor 
of momentum. 

VI. DISCUSSION AND CONCLUSIONS 

Here we do a comparison with the limits obtained in Ref. jfj], in which limits from relic 
density, Tevatron, and gamma-ray are shown. Comparisons are summarized as follows. 

1. The limits due to relic density obtained in this work are consistent with results of 
Ref. [6J. 

2. In Fig. 5 of Ref. [6], the limits for their -Di_4 (corresponding to our O7_io, Dirac DM 
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TABLE II: The 3a lower limits on the operators listed in Table HI We take the coefficient C = 1 
with m x = 50, 100, 200 and 400 GeV. We have used the PAMELA data points above the kinetic 
energy T = 4 GeV in our analysis, because of the large uncertainty of the theoretical background 
at low energy. The x 2 (bkdg) = 5.0. 



Operators 


A (TeV) 






m x (GeV) 


= 50 100 200 400 




Dirac DM, 


Vector Boson Exchange 




Oi = 


= (x7 M x) {qi^q) 




1.15 1.34 1.57 1.66 


2 ~- 


= (X7 M 7 5 x) (<Z7/W) 




0.033 0.038 0.045 0.047 


o 3 -- 


= (X7 M x) (<?7m7 5 <?) 




1.15 1.34 1.57 1.66 


A ~- 


= (X7 M 7 5 X) (97^7 5 9) 




0.19 0.15 0.11 0.09 


o 5 -- 


= (X^X) (q<rp,q) 




1.37 1.60 1.87 1.97 


o 6 = 


= (x^ u i 5 x) (?<v?) 




1.36 1.60 1.87 1.97 




Dirac DM, 


Scalar Boson Exchange 




o 7 -- 


= (xx) (qq) 




0.012 0.013 0.014 0.015 


Os-- 


= (x7 5 x) (qq) 




0.12 0.13 0.14 0.15 


o 9 -- 


= (xx) (qi 5 q) 




0.012 0.013 0.014 0.015 


O10 


= (X7 5 X) (<?7 5 «) 




0.12 0.13 0.14 0.15 


Dirac DM, Gluonic 


On 


= (XX)^,G^ 




0.013 0.015 0.019 0.027 


On 


= (X7 5 X) 




0.13 0.15 0.19 0.27 


O13 


= (xx)G tlv G^ 




0.013 0.015 0.019 0.027 


O u 


= (X7 5 X) G^yG^ 




0.13 0.15 0.19 0.27 


Complex Scalar DM, Vector Boson Exchange 


Ol5 


= (x^mx) (qi^q) 




0.033 0.038 0.045 0.047 


Oib 


= (x^x) (<?7 / V<?) 




0.033 0.038 0.045 0.047 


Complex Scalar DM, Scalar Vector Boson Exchange 


17 


= (x f x) (99) 




0.16 0.13 0.099 0.074 


Ois 


= (x f x) (<?7 5 <?) 




0.16 0.13 0.099 0.074 


Complex Scalar DM, Gluonic 


Ol9 
O20 


= (x t x) G^ V G^ U 
= (X ] X)G^G^ 


14 


0.18 0.15 0.15 0.18 
0.18 0.15 0.15 0.18 



with (pseudo) scalar-boson exchanges) are shown. The limits obtained from FERMI 
gamma-ray are about the same as what we obtained from PAMELA antiproton data. 
The limits from FERMI improve with increasing m x while it is almost flat in our case 
(see Fig. ED. 

3. In Fig. 6 of Ref. js], the limits for their -Ds_s (corresponding to our Oi_4, Dirac DM 
with (axial) vector-boson exchanges) are shown. The limits on obtained from 
FERMI gamma-ray are about 0.1 — 0.5 TeV while the limits that we obtained from 
antiproton data are 1.1 — 1.7 TeV, significantly stronger. The limits on 02,4 are velocity 
suppressed and are 0.03 — 0.2 TeV. 

4. In Fig. 8 of Ref. [6], the limits for their (corresponding to our On^s, complex 
scalar DM with scalar-boson exchanges) are shown. The limits that we obtained from 
antiproton data are slightly stronger than those from FERMI gamma-ray data. 

5. In Fig. 9 of Ref. [6], the limits for their C 3i4 (corresponding to our Oi 5j i 6 , complex 
scalar DM with vector-boson exchanges) are shown. The limits that we obtained from 
antiproton data are weaker than those from FERMI gamma-ray data. 

6. The limits from Tevatron data [6] are rather insensitive to m x , except when m x > 
200 GeV. The limits obtained from FERMI gamma-ray data roughly improve with 
increasing m x |6|. The limits obtained from PAMELA antiproton data are in general 
quite flat. 

In summary, we have used an effective interaction approach to investigate the effects of 
dark matter interactions with light quarks on antiproton flux from the Galactic halo. We 
have assumed a standard halo density and used the GALPROP to calculate the diffusion. 
The obtained antiproton flux is then compared with the PAMELA data. We have success- 
fully used the data to obtain a 3a limits on the scale Aj. The best limits are from the Dirac 
DM with vector-boson or tensor-like exchanges. The limits are about 1 — 2 TeV. While 
the other operators give milder limits. Note that these limits from antiproton flux are lower 
limits on Aj. With the requirement of not exceeding the relic density of the cold dark matter 
deduced from the WMAP, we also obtain the upper limits on Aj. Therefore, both the relic 
density and antiproton constraints give a valid range for each Aj, e.g, 1.6 TeV < A x < 3 TeV 
for m x = 200 GeV. This is a very useful piece of information on the effective interactions of 
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dark matter with the SM light quarks that can give useful implications for collider searches 
and direct detection. 
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Appendix A: Annihilation Cross Sections 



Here we list all the formulas for annihilation cross sections of the operators 0\ to O 
da, 1 N c f3 Q 



[u 2 m + t 2 m + 2s(m 2 x + m 2 q )] , 
[u 2 m + t 2 m + 2s{m 2 q - m 2 x ) - 8m 2 q m 2 x ] , 
[u 2 m + t 2 m + 2s{m 2 x - m 2 q ) - 8m 2 q m 2 x ] , 
[u 2 m + t 2 m - 2s(m 2 x + m 2 ) + lQm 2 q m 2 x ] , 
[2(u 2 m + t 2 J + 2s(m 2 x + m\) + 8m 2 q m 2 x - s 2 ] , 

rn 2 qNc X 



dz ' A 4 16tts/3 x 

da 2 1 N c P q 

dz ~ A 4 16vrs/3 X 

das . 1 Ng P q 

dz ' A 4 16vrs/3 X 

dai 1 Nc Pq 

dz ~~ A 4 16tts/3 x 

da 5 1 Nc P q 

dz ' A 4 4tts/3 x 



dz .V' 32tt ^« 

c 



da 8 m 2 N c sffi 



dz A 6 32vr /3 X ' 

dag _ ^N^n a 

dz A6 32tt t xf g ' 
rfaio N c sP q 



dz ~ A 6 32tt /3 







^ n al l -s 2 P x , G 



cfe A°32tt 

dayi a 2 1 s 2 , 
= A* 32^ ' ^ 

(ut — m 2 (u + t) — m x + m q ) , (, 
(ut - (m x - m 2 q ) 2 ) , (, 





(io"i3 dan 

dz dz 

da 14 da yi 

dz dz ' 

dai 5 1 AT C 

~ A 4 47rs/3 X 

^16 1 N c Pg 

dz ~~ A 4 4irsP x 

dan ™ 2 q N c P 3 Q 

dz A 4 16vr/3 X 

dais . m l N c Pq 

dz A 4 16vr P x 

da±g a 2 1 s 







dz A 4 16vr /3 X ' 
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where s is the square of the center-of-mass energy, z is the cosine of scattering angle, u m = 
u — m 2 x — m 2 , t m — t — m 2 x — m 2 , (3 X — (1 — 4m^/s) 1 / 2 , /3 q — (1 — Am^/s) 1 / 2 , and N c is 
the color factor (3 for quarks). We have set the coefficient C = 1 in these formulas. The 
annihilation rate in the nonrelativistic limit will then be given by o ■ (2(3 X ). 
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